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Abstract. We prove that the singularities of a potential in the two and three dimensional 
Schrodinger equation are the same as the singularities of the Born approximation (Diffraction 
Tomography), obtained from backscattering inverse data, with an accuracy of l/2 — derivative in 
the scale of L 2 -based Sobolev spaces. This improves previous results, see _28_ and [181 . removing 
several constrains on the a priori regularity of the potential. The improvement is based on the 
study of the smoothing properties of the quartic term in the Neumann-Born expansion of the 
scattering amplitude in 3D, together with a Leibniz formula for multiple scattering valid in any 
dimension. 



The inverse scattering problem for Schrodinger potentials deals with the uniqueness, reconstruc- 
tion and stability of the potential q in the Hamiltonian H = A + q from the far field pattern of the 
generalized eigenfunctions or scattering solutions. These are the unique solutions of the asymptotic 
boundary value problem 



where the function u out satisfies the outgoing Sommerfeld radiation condition, which means, for 
compactly supported potential q, that u has asymptotics as |ac| — > oo 



1. Introduction 
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x/\x\. 
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The function A(9' , 9, k), k > 0, 9, 9' in the unit sphere S n 1 , is known as scattering amplitude or 
far field pattern. 

The inverse problem for whole data is formally overdeterminate, as one easily can see by counting 
variables. For this reason, to avoid redundancies, some kinds of partial data are selected for the 
inverse problems. The selection of these data is motivated by numerical experience and applications. 
The most celebrated sets of partial data are the following: 

• Fixed energy data. We assume as data the values A(9' , 9, k) for fixed k and free 9, 9' E S™ -1 . 
Uniqueness of the inverse problem in this case was studied by [16], [17], [23], [32] ■ The 
approach to this problem is related to the Calderon- Sylvester- Uhlmann complex exponential 
solutions, used in the Electrical Impedance Tomography inverse problem. The stability 
happens to be very weak. 

• Fixed angle data. The knowledge of A(9',9,k) for fixed 9, free 9' e § n_1 and free k > 
is assumed. The uniqueness of the inverse problem is open and only generic and local 
uniqueness is proved under a priori regularity assumptions on the potential, see [30j . 

• Backscattering data. One assumes A{— 9, 9, k) for free 9 £ § n_1 and free k > 0. The 
uniqueness of the inverse problem is only proved generically and for small potentials, see [7], 
see also [22], [30], [H]. 

In practical applications the actual potential is substituted by the so called Born approximation 
of the scattering amplitude. The procedures to imaging the Born approximation from the scattering 
data are known as Diffraction Tomography. 

The different Born approximations are obtained, in the frequency domain, from the formula 

q apP rox(0 = A(u>,9,k), 

where £ is given by the redundant relation 

£ = k{u)-Q). 

If 9 is fixed (fixed angle data), we use the change of variable £ = &e(k,uj) — k(u — 9) to define the 
Born approximation = A(ui, 9, k). Notice that this change of variable becomes singular on the 
hyperplane £-9 = 0. 

For backscattering data the Born approximation is given in the frequency domain by the polar 
coordinates 

q B {-2k9) = A{-9,9,k). (1.3) 

This fact makes backscattering data more natural and simpler than fixed angle data for diffraction 
tomography. 

The use of the Born approximation is not, in general, justified on a mathematical basis: one would 
like to know how much information on the actual potential q is contained in the Born approximation. 
This problem has been treated by several authors. In full data case see [19], [21] . [20] and [2]. 
For fixed angle data and backscattering data, both of which are formally well determinate, the 
justification of diffraction tomography was studied in [26] (fixed angle) and [11], [18], [28], [24] 
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(backscattering) . We would like to remark that each of the last two types of data require the 
analysis of special multilinear operators which are not related. 

In this work we study the case of backscattering data in dimension two and three, we continue 
and complete the research of [IB], [Hj and [H], by removing some constrains in their results. We 
prove that the diffraction tomography is a migration scheme, see [3], within an accuracy of at least 
1/2. This is to say that the most singular parts of the actual potential can be reconstructed from 
the Born approximation up to a certain order (the accuracy of the migration) . The determination 
of this accuracy is very important to design numerical methods, adapted to the spaces in which one 
expects to obtain the information on the actual potential from real scattering data. We prove 

Theorem 1. Assume that n £ {2,3}, q is a compactly supported function in W a ' 2 (M. n ) and a > 0. 
Then q-q B G W /3 ' 2 (R n ) + C°° (M™) , for any P £ R such that < < a + \ . 

In Theorem [1] the regularity is measured in the scale of L 2 -based Sobolev spaces. The optimality 
of this accuracy in this scale of spaces is, so far, an open and interesting question. 

The procedure to justify the migration scheme is to study the smoothing properties of the multi- 
linear terms in the Neumann-Born expansion of the scattering amplitude (multiple scattering). 

Physical evidence suggests that multiple scattering is strong in the case of backscattering data. 
The control of double and triple scattering in 3D, within an accuracy of 1/2, was obtained in [2"8] . 
but their result together with the general estimates for multiple scattering do not suffice to assure 
that, for a potential q a priori in the Sobolev space W a ' 2 , the error q — qs is in W 13 ' 2 for any 
(3 < a + 1/2; the restriction < a < 3/4 is needed. In the range a > 3/4 known estimates 
of quadruple scattering became worse than those of double or triple scattering. The study of the 
accuracy of the Born approximation requires, then, to improve the estimates of the quartic term in 
the series. We accomplish this in the present work. We also extend the results, which previously 
were only studied for a < 3/2 in 3D and for a < 1 in 2D, to any a > by using a Leibniz' type 
formula for the derivatives of multiple scattering terms (see §C.l in [25 ). 

In dimension three, we only are able to prove that the errors due to double, triple and quadruple 
scattering are a half of a derivative better than the actual potential, as opposite to the 2D case 
where the regularity increases with the order. 

Result from [28] , [24] together with Corollary Q] allow us to state the following result concerning 
reconstruction of classical discontinuities from backscattering in 2D: 

Theorem 2. Let q compactly supported in W a,2 (M. 2 ), where a > 0. Then q — qs is a continuous 
function. 

In fact, it was proved (Theorem 2 in 28J), that, for such a q, the quadratic term is a continuous 
function. Holder continuity of the cubic term is obtained since it is in W 13,2 for all /3 < a + 1, see 
Theorem 1 in [21]. The remainder is controled by Corollary [1] 

In the three dimensional case, it follows from Theorem 1 that the whole non continuous part of 
the actual potential can be reconstructed from the Born approximation, assuming a priori that q is 
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in the Sobolev space W a ' 2 for some a > 1. Notice q might have some discontinuities if a is between 

1 and the 3D Sobolev exponent 3/2: 

Theorem 3. Let q compactly supported in W a,2 (M. 3 ), where a > 1. Then q — qs is a continuous 
function. 

From the previous work [28J it follows that in 3D the discontinuities in the case of a piecewise 
regular potential can be reconstructed from the Born approximation (the result is not stated in 
[28] but it is similar to Corollary 0.1 in [18] in the 2D case). By using the evolution equation the 
reconstruction of conormal singularities was achived in [11] , On one hand Theorem 3, as far as we 
know, is the first result of reconstruction of discontinuities in 3D, without assuming special structure 
of the singular set but, on the other hand, one expects that q G W a ' 2 , for any a > 1/2 suffices for 
the reconstruction of discontinuities. So far, this improvement has not been achieved. We know 
from Corollary [1] that the high frequency Neumann-Born series for j > 5 converges to a Holder 
continuous function for a > \. 

An important feature of Theorem [1] is the fact that, regardless of the a priori regularity assump- 
tions on the potential, the accuracy of the migration scheme is at least 1/2. This independency is 
important to construct any recurrence scheme, in order to obtain further information on the actual 
potential from scattering data. In the case of fixed angle data, one can define a modified Born ap- 
proximation by inserting the error q — qs in the quadratic form, see [26] . This increases the known 
accuracy for rough potentials q, but an inconvenient to iterate the procedure is the dependency on 
a of the accuracy. 

Finally, we remark that in the higher dimensional case, the order of accuracy is an open question. 
We believe that 1/2 also applies, but the technical complexity of our approach makes it necessary 
to look for a new point of view on the problem. The treatment of the 3D problem due to Lagergren 
[14] . [15] , based upon a time dependent expansion of the backscattering operator, also requires a 
very technical treatment of its multilinear term. See also [I] and [5]. 

Notation and definitions. We will write J-f or / to denote the Fourier transform of /. T~~ x de- 
notes the inverse Fourier transform. The letter M denotes the Hardy-Littlewood maximal operator. 
We denote the two-dimensional Hausdorff measure in R 3 by a. The expression |a;| ~ 2 _j '|jj| refers to 

2 J 1 1 ^7 1 < M < 2 -j ' +1 |t7|, for j G Z, x,r] € R 3 . We will use the homogeneous and non homogeneous 
Hilbertian Sobolev spaces. With a 6 R, we denote 



W a ' 2 (R n ) : 
W a ' 2 (W l ) : 



{/ G S'(R n ) : | ■ \ a Tf(-) G L 2 (R")}, 

{/ G S'(R") : (1 + | • | 2 ) Q/2 ^/(-) G L 2 (R")}. 



Let <q, £ G R 3 \ {0} . We write 




(1.4) 



refering to the sphere centered at | and radius ^ and 

A(0 : {■>■■ O 0} 



(1.5) 
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denotes the plane orthogonal to £ that contains the point £. We denote r(?7) 3 := T(r/) x T(rj) x T(r)) 
and r(?7) 2 := T(r}) x T(r)). 

Let F(rj) given by the integral on a manifold A(rj) of some function. Since our proofs are based 
upon a decomposition of A(rf) in several subdomains D(r)) C A(rj), we will denote by Fjj(rj) the 
same expression when we restrict the integration to the subdomain D(ri). 

The outgoing resolvent operator for the Laplacian is defined, in terms of the Fourier transform, 

by 

R k (fm = Ht\ 2 +k 2 +io)-if(o. 

We define the operator Qj in the following way 

QM(0 = J e ik9 -v(qR k y-\q(-)e ikS -U)(y)dy, (1.6) 

where k = |£|/2, 9 = — With these expressions for k and 9, we define the multi-linear form 
Qj (/i , ... , fj) in the FT side as 

^(QiCfi s .--./i))(0 := / e ike -y{hR k ){hR k ) ... (f^R^f^e^Mdy. 
We denote the high frequency version 

Qi(?):=^V(|- \mOM(-)), (1.7) 

where Qj(q) is defined in JTJJ) and x* € C°°(M) with **(t) = 1 if * > 2C , x* (*) = if t < C , for 
a certain constant Co > to be chosen (see section 4). Notice that the cutoff near the origin allows 
us to reduce the estimates of Sobolev norms to the estimation of homogeneous Sobolev norms. 

We also write Qj(fi, ■ ■ - ,fj) = (x*3~ (Qi(/i; ■ ■ • > fj)) ) ■ We will admit the abuse of notation 
Qj{q) = Qj(q, ■ ■ ■ , q) and Qj(q) = Qj(q, ... ,q). 

The permutation group of order k is denoted by S k ■ For multi- indexes (3 and 7 in N™ , we use the 
standard definitions of j3\, \j3\ and j3 < 7. 

We use the letter C to denote any constant that can be explicitly computed in terms of known 
quantities. The exact value denoted by C may change from line to line in a given computation. 

2. Preliminaries and results 
We obtain the so called Lippmann-Schwinger integral equation by applying the outgoing resolvent 

to CEO) 

u(x, 9, k) = e ikx - 9 + R k (q(-)u{- . 9, k))(x) . (2.1) 
The key operator in the above integral equation is 

T k (f)(x)=R k (q(-)f(-))(x). 

There are several a priori estimates for R k that allow to prove existence and uniqueness of Lippmann- 
Schwinger integral equation. Usually, Fredholm theory applies and everything follows from compact- 
ness arguments, Rellich uniqueness theorem and unique continuation principles, in the case of real 
valued potentials. The solution can be obtained in several situations (these cases do not require q to 
be real) by perturbation arguments, assuming that the energy is sufficiently large, k > fcg > , where 
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k depends on some a priori bound of the potential q . As an example we may consider compactly 
supported q G L r (R") for some r > § . In this case, which is the one considered in this work, the 
resolvent operator Rk is bounded from L p (M. n ) to LP (M n ) with norm decaying to as k — > oo when 
i — y = ~ , see pQ , [13] and see also [25] . This together with Holder inequality proves that for big k 
the operator Tfc is a contraction in LP and then existence and uniqueness of solution of (|2.1[) easily 
follow and u can be expressed as a convergent Neumann-Born series. 

Once the scattering solution is obtained we may prove that the far field pattern can be expressed 

as 

A(6',6,h)= [ e- M '-yq{y)u{y,e,k)dy, (2.2) 
see [7] where this is used as a definition for non compactly supported potentials. 

By inserting the series u in (|2.2[) one obtains the Neumann-Born series of the scattering amplitude 
for k large enough (high frequency Born series): 

oo 

A(0',6, k) X * (k) = q(k(0' - 0)) X * (k) + ]T qj(q)(e', 6, k) , (2.3) 

3=2 

where 

q j (qW,6,k) = X*(k) f e- lkd ' ' ' y {qR^- 1 {q{-)e lkd ^){y)dy , 
and x* is a cutoff function near the origen (see the notations). 

We deal with the backscattering inverse problem, for which one assumes the data with the direc- 
tion of the receiver opposed to the source direction (echoes), i.e. A(— 6, 6, k) . The inverse problem 
is then formally well determined. In this case the Neumann-Born series for the scattering amplitude 
is 

oo 

A(-e, e, k) x * (k) = q(Ox* (k) + Qi . ( 2 - 4 ) 

where £ = —2k0 and the j-adic term in the Neumann-Born series Qj{q) is given by the operator 
(|1.7[) . We define the Born approximation for high frequency backscattering data as 

q^5(Z) = A(-6,0,k) X *(k) 

where £ = —2k9. Notice that the series (|2.4[) is addapted to the reconstruction of singularities, since 
<Zb, h — qs and q — J-^ 1 (q(-)x* {\ • |/2)) are C°° functions. 

We denote the remainder term in the high frequency series as 

oo 

R-K<7) = E^(<?)- 

3=1 

The main part of this work, which is §3, is due to obtain the control of the term Qi(q) in dimension 
three: 

Theorem 4. Let us assume that q is a compactly supported function in W a ' 2 (M 3 ) , for < a < 3/2 . 
Then Q 4 (q) € W@ ' 2 (R 3 ) , for any f3 such that < f3 < a + 1/2 . 

We also prove in section §4: 
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Theorem 5. Let n G {2,3}, q € W a ' 2 (M. n ) compactly supported and < a < n/2. Assume that 
Co > max{(2|jg||i4/o,2) , 1} , I = 4 if n = 2 and I = 5 if n = 3, see |1.7| j. Then, for any (3 € K 
smc/i i/iaf /3 < cf + 1 £/ie remainder in the high frequency Born series R; converges to a function in 
W^< 2 (R"). 

From Sobolev embedding theorem we obtain, 

Corollary 1. In the hypothesis of Theorem [5l assume also a > in 2D and a > | in 3D. Then R; 
is a Holder continuous function. 

Theorem 1 in the case < a < n/2 will follow from the above theorems, together with the 
following estimates for the quadratic and cubic terms, see [28] and [2"4"] : 

\\Q2(q)\\w^<C\\q\\ w -iA<l\\w^, (2.5) 

WQMWw^ < c + U\lAW\\ w -^ + h\\ w -i,h\\w-^) h\\w^ > ( 2 - 6 ) 

where the dimension is n = 3, (3<a+l/2,e:=a + ^ — /3 > 0. For dimension n = 2we have (|2.5p 
for /? < a + 1/2 and 

||Qa(g)||*w < C (||9|U a ||g||^_ 4 . a + ||g||! 2 ) \\q\\ wa ,2 , (2.7) 

when < < a+ 1. 

In §5 we give the procedure to extend the above results to the case a > n/2. The key is Theorem 
[6] which is a Leibniz' type formula for derivatives of multiple scattering terms. 

The proofs are very involved and technical. For this reason, we only include the details of the 
proof in the key case of Theorem [H see Proposition [3] in §3.1. In other cases, we just sketch the 
proof and try to convince the reader that similar arguments work. 

3. Proof of Theorem [H 

Theorem U follows from the estimate 
Proposition 1. For q,a under the same hypothesis of Theorem^it holds 

WQMWw^ < c (h\\h + Ikll^alMli* + h\\w-h*Mw-^hh* (3.1) 
+ IMIi»IMItf,-i-.. a ) h\\w^, (3.2) 

for all /3 68 such that < < a + 1/2 , where e:=a+|-/3>0. 

The quartic term in the Neumann-Born series for backscattering data is given by 
Q^qj(0~ [ e lke -y(qR + {k 2 )Y{q{-y ke ^)(y)dy : 

JR 3 

for any (6l 3 , where £ = — 2k6 , that is, k = ^ and — — t|t . From Lemma 3.1 in [27] , this term 
can written as 
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Proposition 2. For any dimension n and r\ € R ra \ {0} , 

rfTv \ f f f q(0q(v - t )q( t - 4>)q{<t> - , c , , q o\ 

QM(V) = P-V. I I I [ M „-0][r-(,-r)] ^ 

+ 2^ P ,.j f [ j( y ( r T) tm P^ MOdrd^ 0.4) 
|??| V y R n J r(n) [r • [rj - t)\ [</} ■ (77 - </>)] 

+ T~\ p - V - / / m cvTT — 1 \i d^dTda(cl>) (3.5) 

-2-—^p.v. / / — d£,do{T)do(4>) (3.6) 

\V\ Jr(r,) Jr( v ) M £ • W - 



Xj?) <M?7 - 0) 

9(0«(»7 - ^ - ~ <M£)*r(T)dcr(0) . (3.. 



\V\ Jr- Mv) Mn) V'iv-V) 



M 3 Jr( v ) Jr(ri) Jr( v ) 

The key to understand the structure of the quartic term is the pure spherical measures part (|3 .8[) . 
Hence we define, for any 77 G R 3 \ {0} 

Notation: 

Q(q)(v) ■= ,4 / / / mq(v-T)q(T-m(<f>-0M0do-(T)do-(<f>). (3.9) 

l 7 7l 7r(n) -'r( J) ) Jr( v ) 

We prove in section §3.1: 

Proposition 3. Let q <G W a ' 2 (M. 3 ) be a compactly supported function with < a < 3/2. Then for 
all (3 < a + 1/2, 

\\Q{q)\\w^ < C (\\q\\h + \\q\\ w -±,2\\q\\h + \\q\\ w -^Ai\\w-^\\q\\^ 
+ \\q\\h\\q\\w-i-c^) h\\w«.» • 

where e := a + ^ — /3 > and the constant C > just depends of a, f3, e and i/ie support of q. 

Now we sketch the estimates of principal value terms (|3.3|) - (|3.7|) We use a decomposition of the 
space into diadic shelves, as it was done for the cubic term in 2D, see [24], and for the quadratic 
and cubic terms in 3D in [28]. More detail can be seen in [25] . 

Let us state, as a model, the main features to control the principal value term Q'(q), given by 
HQ'(q))(v):=^ P .v.f f I m ^-^-^-^ da(0da(rW. (3.10) 

M Jw Mr,) Mi) 4>-{ r n-4>) 



The key to estimate this principal value operator is to control the term: 

1 

M 4 Jr s ( v ) Jr( v ) Jt{ v ) 



F(Qs(q))(ri) ~X(i-Soo)(l»?l)o4 / . / /_. J(0q(v-r)q(r-<f>)q(4>-0 da^)da{r)dct>, (3.11) 



where 

r s (v) := € M 3 : ||0 - v/2\ - |r?|/2| < ^lr?|> (3.12) 
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Comparing (|3.11|) with (|3.9p . we observe that we replace the sphere T(rj), in which the variable 4> 
runs, by its tubular neighborhood Ts{rj) of width 5\r)\ in the normal direction. Notice that da n {<j)) = 
lira ^j— - Xr s (ri) {4 1 ) d(f>, where da v {4>) denotes the measure on the sphere T(rj) induced by Lebesgue 

measure d<f>. For 5 small, we have that Qs{q){Tj) ~ X( < 5- 1 .oo)(l 7 ?l) ( ^Q('?)( 7 7)- I n this way we may 
expect estimates for the Sobolev norm of Qs(q) obtained from estimates of Q(q) multiplied by S. If 
one follows the lines of the proof of Proposition [3J one gets the following 

Lemma 3.1. Assume that q £ W a ' 2 (R 3 ) is compactly supported and < a < 3/2. Let (3 < a + 1/2 
and e ~ a + 1/2 — /3 > 0. Then there exist S± > and 7 = 7(e) > 1 such that for all S satisfying 
< S < 5± it holds 

\\Qs{q)\\w^ < (|M|| 2 + \\q\\ w - h A\l\\h + U\ w - h M\w-M^ (3-13) 

+ lkll! 2 |kll^-^e, 2 )l|g||^, (3.i4) 

where C just depends on a, f3, Si and the support of q. 

Now, to estimate the term (|3.10[) . we use a decomposition of the Euclidean space R 3 in a similar 
way as was done in [24] for 2D: 

P°g 2 \V\] 

R 3 = r ir (77)u |J rjWu^W. ( 3 -i5) 
3 =31 

where j\ is the lowest integer such that j\ > 1 — log 2 ((5i) with <5i from Lemma f3. It |?7| > 2 J1 ~ 1 and 
I\- (rj) :—{(/) e R 3 : \ \(f> ~ r,/2\ - \r,\/2\ > 2^ +1 |ry|}, 
Tjir,) := {cj> G R 3 : ||0 - V /2\ - \ V \/2\ ~ 2^\ V \}, j x < j < [log 2 \ v \], 
I^fo) := {<f> G R 3 : ||0 - 77/2| - \n\/2\ < ^M}. 

Remark. Technically this partition only makes sense for ji > 3, but this is not a constraint if we 
demand Si < 1/4, since ji>l — log 2 (5i). Notice that 1^(77) C Too (77), where 

^ooiv) :^{cf>eR 3 :\\cb- V /2\-\ V \/2\<l}. 

This decomposition is used to split the operator (|3.10p . To control the operator corresponding to 
the annulus terms Lemma l3~Tl with S = 2~ J+1 , suffices. To deal with the central term, corresponding 
to (77), which is close to the singularity T(r]), we use again Lemma 13. II and the following 

Lemma 3.2. Let fj € W a,2 (R 3 ), j = 1, . ..,4, be functions such that fi, / 2 are compactly supported 
and < a < 3/2. We denote 

HQUhJ2j3,h))(v) -=X*(V)A [ I I \M)h(v~r)\ (3.16) 

x \f 3 (r - 4>)h^ ~ 0\MOMr)d<t>. (3.17) 

Then for any f3 < a + 1/2, 

1 1 ^ oo 

< C(a,/3,SUpp/i,SUpp/ 2 ) ( ^2 ||/<r(l)|U» \\f<r(2) |UHI/<7(3)|U=|l/<7(4) \\w<.,2 ( 3 - 18 ) 
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+ E ll/r(l)ll^-i, 2 |l/r(2)|UHI/r(3)|UHI/r(4)Hw^ (3-19) 



T6S4 



2J IIA»(i)llw-i. 2 11/^(2) llxv-=,2 11/^(3) Hi 2 11/^(4) II vv-. 2 (3.20) 



<w 



E ll/p(l)llw-i-^2||/p(2)llL 2 ll/p(3)lli 2 ll/p(4)llvy°.2 ) 1 ( 3 - 21 ) 



P6S4 



where e := a + ^ — /3 > 0. 



To estimate this central term, dealing with the principal value, one needs to use the cancelation. 
We must replace the integral on the ring (77) by J*r/+(jj) + /r~(rj)' wnere 

r+(T?) = {e e m 3 : £ < ie - vm - h/2 < 2-^ '"'^mi 

r-(^) = {( £ I 3 : e < h|/2 - £ - V2| < 2-P°fc l""" 1 ^!}. 

The map F : r^~(?7) — > T~ (77), given by symmetry with respect to T(rj) allows us to pass to the limit 
when e — > + . To cancel the singularities we use an estimate, due to Calderon, for first differences 
in terms of the Hardy-Littlewood maximal operator M (as in [24], several standard reductions are 
also needed): 

Lemma 3.3 (see [12]). Let u £ W 1>p (SL n ), p > 1, a G E™. Then 

\u{x) -u{x-a)\<C \a\ [M{Vu)(x) + M{Vu)(x - a)}. (3.22) 

After some changes of variables in the integrals involving F, we reduce to study the following 
terms: 

1 



r , , , \q(Oq(v-T)q(T-<t>)\MVq(<t>-Qda(Oda(T)ckt>, 

\V\ JroaWJrwJrfr) 

r4 / II \mQ{v-T)\Mq{T-4>)MVq{ct>-t)da{£,)<kT{T)d<t>, 
lis / / / \m^n-r)\MVq{r-4>)Mq{4>-OMOdcr{T)d4>, 
Ha / / / l9(O«(»7-r)|MVg(r-0)|g(0-OI^(Odff(T)d0- 

The proof of Lemma 13.21 follows the lines of the proof of Proposition [3l Heuristically, Lemma 
3.21 is derived from Proposition [3] replacing the domain T(r]) for the variable <f> by the tubular 
neighborhood (rj) which is the result of widening the sphere T(?7) a distance 1 in the normal 
direction. Nevertheless there is an additional difficulty which has to be managed: the fact that 
neither f^, nor f± are compactly supported and their Fourier transform can not be controlled by the 
maximal operator using Lemma 16.21 But we must keep in mind that we can apply Lemma 16.21 to 
two functions, /1, /2, which are compactly supported, and the integral of I/3I 2 or I/4I 2 in (f> can be 
bounded by the L 2 -norm using that the variable <fi is solid. After these comments, we omit the long 
and tedious proof of Lemma 13.21 The reader can see all the details in a similar situation for the 
cubic term in 2D (Lemma 2.2.3 of [2"5]). 



The key to control the principal value term (|3.3|) remains in the following lemma whose proof 
follows the lines of the proof of Proposition [3] 
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Lemma 3.4. We denote 



x / / \mttv-TMT-4>)q(<t>-t)\d£dTd<l>. 

Jr Sl ( v ) Jr S2 ( v ) JT l3 ( v ) 

Let q and a as in Lemma \3.1\ and (3 < a + 1/2, e = a + 1/2 — [3 > 0. Then there exist 8q > and 
7* = 7*(e) > 1 such that for all 61,82,63 satisfying < 81,62,63 < 80 it holds 

\\Qs u 5 a ,5 s (q)\\wf>.* <C(8i8 2 83r \\q\\w^^), 
where C only depends on a, [3, 80 and the support of q. 

Analogously to the comment about Lemma \'6. II above, this result should not be surprising since 
FQs 1 ,5 2 ,8 3 {q){ri) ~ 8162833- 'Q(q)(f]), for 81,82,83 small. To estimate the term (J3T3J) we have to take 
the partition (13.15|) of M 3 with ji the lowest integer such that j\ > 1 — log 2 (<5o)> for 8q from Lemma 
13.41 In particular, the control of the ring terms 

/ / / TFT FTT — 1 \ul~i -TT^dTdcf), ji < J,k,l < log 2 77 

follows from Lemma l3~4l with Si = 2~' J+l ,S 2 = 2~ k+1 ,8 3 = 2~ l+1 , together with the fact that 

\C-(v-0\ = (\C-l\ + \l\M-l\-\l\)>c\r 1 \ 2 8 1 , 
where c > and we use definition (|3.12[) . 

To estimate the central term 



'r^W Jr^v) JrZeW (»? ~ ' (v ~ r)][<f> ■ (r? - 0)] 
we must replace each integral on the ring 1^,(77) by J r +/, ) \ + Jr~(»j)' then, use the ma P -F : 

r+(?7) — » r~(ry). We need again Calderon estimate for first differences and its analogous estimate 
for second differences: 

Lemma 3.5. Let u G W 2 ' p (R n ), p > 1, a, b, c G K". TTien 

4 

+ 6) -u(x) -u{x + b-a)\ < C\a\ \b\J2 M 2 \D 2 u\(xj), 

,?'=i 

where D 2 u denotes the matrix of derivatives of order two and x\ = x, x 2 — x — a, x% = x + b, 
X4, = x + b — a. 

These tools allow us to reduce to a sum of integrals, analogous to those written after Lemma 3.3 
for the case p.7p . 

3.1. Proof of Proposition [3j 

Let us split the set r(?y) 3 into the following regions 

Or,) :={(£,T»Gr(r7) 3 : |*-£|>^, 1^1 > 
H(V) :={(e,r,0)Gr(^) a : \cf> - £| > , |0 _ T | < J1J 
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nm :={K-^)er(,f : |0-£|< Jl,|0- r |>ili 

IV(V) :=ht,T,<fi)er(r,)* : |0-f|<Mj _ T |< JjjLJ. 
In this way, we can write Q(q) = Qi(q) + Qn(q) + Qin{q) + Qiv(q)- We will prove that 

\\Qi(q)\\w^ < C + IMI^-^lklli") h\\w»-v*.*i ( 3 - 23 ) 

\\Qn(q)\\w^ < C + H\ w -^h\\h) hWfre-w,*, (3-24) 

IIO/v(g)||^.. < c (lkll^-i, 2 |klU-e, 2 |klU= + ||g||! 2 |kii^_i_ £ , 2 ) |M|^_i +e , 2 , (3.25) 

provided that e > 0. Note that Qin(q) satisfies the estimate (|3.24[) since Qn(q) = Qni(q)- 

, we have 

\mmv - r)q(r - 0')g(0' - 0\d*(0da(T)da(<t>') 



Proof of estimate (|3.24p . Taking the change of variable = r) ■ 
1 



Qn(q)(v) = tti3 



n( v ) 



1 

" w 

We decompose 



{(t,T,<t,):(tT,T,-<t>)£II(r,)} 



W(0<1(V ~ r)q(r + - 77)9(77 - - o\ M0Mt)M<I>) ■ 



{(e,r»er(^) 3 : ^,T, V -4>)eII{ V )}=\J(II k ( V )UlI k (r 1 )), 

fc=l 

where for any k e N, we denote 

:= ||r? - r - 0| < K , \ v - £ - 0| > K , |0 _ £| ~ 2 - fc h| , |0| < Id 

:= {to - r - 0| < igL , _ f _ 0| > K , |0 _ f | „ 2 - fe h| , Id < |0| 
with (£, r, 0) e F(r}) 3 . We have 

and then to prove (|3.24p we use 

IIQ//(<7)llw//3,2 < E (ll QiiMWwp.z + \\Qji k (<i)\\wf>.*) ■ 



k=l 



+00 



k=l 



For each fe > 1 we claim 

WQnMWw^ <C2- k i 2 \\qf L2 \\ q \\ wfi _^ 2 , 

\\Q rik (q)\\ w ^<C2- k ^ (\\q\\ L 2 + \\q\\^_ h2 ) \\q\\hh\\^-^ 



(3.24a) 
(3.24b) 



In the following, we use the notation in Lemma 16.31 which is the key of the proof of the above 
claims. 



Proof of claim (|3.24ap . By Cauchy-Schwartz inequality, 

|g(O9(T + 0-^)| 2 ^)d<x(r)d ( 7(0) 



n k (n) 



(3.26) 
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\q(v - r')q(v - 0' - OfdaiOda^daW) . (3.27) 
iikin) ) 

If we widen the sphere V(r\) until ri(?7) :— (x G R 3 : |a; — ^ | — ^ < l| , by part (1) of Lemma 
we have 

T + <j)-r!)\ 2 d(j{T)da{(j))<C ( [ Mq{x + cj)-rj) 2 dxda{(j)) 

where the last inequality follows from the boundedness of Hardy-Littlewood maximal operator in 
L 2 (M 3 ) and Plancherel identity, since the measure of T(rj) is | ?7 1 2 - In the same way, 

\q( V -r')\ 2 da(T')<C\\q\\ 2 L2 , (3.28) 

r(v) 

If (£, r, (ft) 6 Hk(v) then |£| > 2~ fe ~ 2 |?7| , and changing the order of integration in £ and rj by Lemma 

urn it holds 

IIQi/ fc (?)ll^. a <C||«||t a / M 2/ W |<z(£)| 2 

\q{v - <t>' - O\ 2 da(0da{4>')da(^ 

A k (n) 

c\\ q \\ 4 L J \m\ 2 [ ^m 2/3 - 4 

A* fa) 

fe II _,|| 4 / \ t\\^ J? ( C\ At r i~ ^ II «ll 6 



<C2%|| 4 L2 / |g(0r^(0^<C2- fc ||5||^||g|K,_ ii2 , (3.29) 
where the last inequality follows from part (i) of Lemma 15731 and fft(£) is defined in (|6.ip . Also, 

^(•?):={(f^)er(f/)xr(„) : |£'-0'| < 2- fc+1 | 7? |, |r?-^-n > i^} . (3.30) 



□ 



Proof of claim (|3.24b[> . We take II k (v) = n l(v) U , where 



II 

2 

II 



\{r,) :={(£,T,<f>)eII k (ri) : \v - 4> - r\ < 2- k - 3 \ V \ } , 
^(77) :={(£,t,0€ JI*(»j) : \ V - 4> - r\ > 2~ fc - 3 M } . 



Let us start with the domain 

-l 



IIM = {(£, r, 0) € L(77) 3 : |t? - r| < 2- fc - 3 |r,|, \ V - £ - r| > |r?|/100, |<£- £| ~ 2- fe |r,|, |£| < |0|}. 
On this region |?y — r| > 2~ fe_3 |?7| holds. By Cauchy-Schwartz inequality, 



\v\" 



Qn] (?)(»?) < ri ( / / } WOSfa - r)| 2 d ( x(e)da(r)da(0) ) (3.31) 

7 / + ^-'/)?(iI-^-r)WWW)] ■ (3-32) 
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By part (1) of Lemma [6.21 we have J r ^ l < 7(£)| 2 <^ <J (£) — Clkllis > an d by this lemma and Fubini's 
theorem, for each £' ,4>' € T(tj) in (|3^2|) . 

\q{r' + c^' - V )\ 2 da{r')<C\\q\\l 2 . 

r(v) 

Taking the change £ = rj — r , and changing the order of integration in £ and rj by Lemma 16. 1[ we 
may write 

\\Qr T i(<i)\\ 2 w^< c h\\h I H 2/3 " 4 / l?(C)l 2 ^(C) 



'{Cer(i,) : |c|>2-*-3|i,|} 



Jr 3 JlveA(C): ICI>2- fe - 3 |j?|> Kl 



A k (.v) 

{^eA(c): ICI>2- fc - 3 M} '' 

<c\\ q \\i 2 2 k [ \m 2 m)dt <c2-%\\u q \\ 2 i(2 , 

where the last inequality follows from part (i) of Lemma 1631 and A k (rj) , F&(£) are defined in (|3.30|) . 

(EH). 

2 2 ~2 2 

We go on with the region II k (r/) . Let us split it as follows: II k {rf) — Hk.ai 1 !) U Hkbvu > where 

— -2 r — 2 

JJ 



^.6 



„(»/) :={(f,r^)G/Z fc (»7) : h-r|< |0|} 
(77) := far,® eTll(r,) : |»j-r|>|0|}. 



-2 



On the region II k a (j]) , we know that if |£| > 2 we can follow the lines of the case Hkiv) ■ So, 

~2 ~2 2 

splitting once more as II kia {rf) = U k>a>1 {rf) U i7 fe)0i2 (r?) , where 

Tll^M ■■= {(£,t,0) e/Ifc, (r?) : |£| > 2^1^} 

= {(£, r, 0) e r(77) 3 : 2- fc ~ 3 M < |r? - 4> ~ t\ < |r?!/100, |r? - f - r| > h|/100, 
|^-e|~2- fe H,2- fc H<|e|<|^|, |»?-r|<|0|} 

^...aW :={(&^)G //*,-(»?) : ICI < 2 ^ | ^ | } 

- {(£, r, 0) S r(r,) 3 : 2" fc - 3 M < \r, - 4> - r| < |»y|/X00, |»7 - f - r| > |ry|/100, 
10 " CI ^ 2- fc |r,|, |e| < 2~ fc H, |f| < \ v - r\ < |0|}. 

We may write 

\\Q Tl l a {q)\\w^ < C2- k / 2 \\q\\l 2 \\ q \\ w0 _ h2 . (3.33) 

~2 _, 

In this way, we reduce to the case II k a 2 (?j) , where \cj>\ < 3 • 2 |?y| holds. 
Remark. From the proof of claim 6 in [55] one deduces to the following estimates: 

\q(x-y)\ 2 da(y)da(x)<C2- k \r 1 \ 2 \\q\f » , (3.34) 

r(r ? )2n{| a =|<|y|<2-^ + i|r ? |,|^-vl<4gl)} 

|g(x - y)| 2 da(t/)rfa(x) < C |r?| || ? ||£ a . (3.35) 

r( J ,)2n{|x|<|y|<2-'= + i|, ) | : 2-fc-3|^|<| a; _ a |<M} 
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By Cauchy-Schwartz inequality as in (|3.26jl . (|3.27|) . the fact that Lemma [6.21 implies (|3.28p , ap- 
plying the change ( = rj — t , since 

|^-C)l 2 ^(0^)<CH||g||i a 

rW 2 n{|C|<|0|<3-2-'=|r,|,2-^3|^|<|C- ( ;i|<i2i} 

by (|3.35p , and changing the order of integration in £ and rj by Lemma 16.11 we have 

\\Q rf {q)\\ 2 w ,,<C\\ q \\\J M 2/ W \m?MZ) (3-36) 
"■°- 2 Jw Jr( v ) 

\q(v-<t>'-Z')\ 2 MZ')M<i>')dV (3.37) 



Ak{v) 

<C\\q\\\, [ l -^fF k (OdC (3.38) 
where the last inequality follows from part (z) of Lemma 16.31 and Ak(rj) , Fk(£) are defined in (|3.30[) . 

& 

Let us go on with the domain 

ThM = {(£> r, 0) G r(r,) 3 : 2" fc - 3 M < \ v - <P - r| < |r?|/100, 

h - 4 - r| > h|/100, |0 - 41 ~ 2- fc |r?|, |£| < |0| < ^ - r|}. 

2 -h-£ 1 

So, if (£,r, </>) € II kb {rf) then |^ — r| > 2 fe 4 |r?| . Following the lines for the case II k (r/) , one 
obtains 

\\Qril b (l)Ww^ < C2- k / 2 }\q\\l4q\\^_ h2 . (3.39) 
We conclude the estimate 

\\Qjf k (q)\\w^ < C2- k/2 (||«|U» + Nl^-i, 2 ) IMI^NI^-f , 2 , (3-40) 
and the claim (|3.24b[) . 

This ends the proof of estimate (|3.24|) . 



□ 



Proof of estimate (|3.23p . Taking the change of variable <j> — rj — <fi' , we have 



1 



Q/(9)fa) = 03 / / / \mttv-T)4(T-4/)qW-t)\fr($da(T)M<l>') 



1 



{(«,T,77-^)e/(»7)} 



For 77 G K 3 fixed, we take the decomposition 

(£, r, 0) G r(r,) 3 : \ V - <j> - £| > Jl, |r/ - J> - r| > j^j = h( V ) U U 7 2 (r?) , 
where 

Afa) :={l0-r|<]^, H<|77-r|, |,y - r - 0| > , |ry - ^ - 0| > ^} , 
:={l0-r|<^, W>\v-t\, \ V - t - c/>\ > , \ V - £ - <f>\ > ^} , 
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h(v) ■= (\4>-t\ > — , In — 7- — 01 > — , In — f — 01 > — 
y " \' '- 400'" ^'-100'" 1 - 100 

with (£,T,(f>) e T(t]) 3 . It holds 

QMin) = QMiv) + QtM)(v) + QiM){v) , 

and also, 

\\QM\\wb,2 < WQiMWwfi,' + \\Q Tl (q)\\w^ + WQiMWw^ ■ 
We claim the following: 

IIQ/i^ll^^ciklli'lkll^-i.". ( 3 - 23a ) 

\\QTM\we.><C\\q\\l4q\\ w0 _i i2 , (3.23b) 
WQiMWwo.' ^ C\\q\\^_ h2 \\qf L2 \\q\\^_ h2 . (3.23c) 



The estimate (13.23)) follows from these three claims. In their proofs we use the notation introduced 
in the key Lemma 16.31 located in the appendix. 



□ 



Proof of claim (|3.23ap . On this region we have \rj — r| > ^ . Applying the Cauchy- 
Schwartz inequality as in (|3.31|) - (|3.32[> . since for j] e M 3 , 4>' e T(r]) fixed, by Lemma ET2l it holds 
f m MOfdaiO < C \\ q \\\ 2l f m \q( V d')\ 2 d<r(C) < C \\qf L , and a(T(r,)) = ir\ V \ 2 , we obtain 

&M)(v)<^h\\h( f \q(T,-r)\ 2 da(r) 
\V\ Vy{Ter(,):|^-T|>i2i} 

W + 4>' ~ri)\ 2 da(r')da{<t>') 

A( V ) 



where 



.1(„): Kr'i'lefM 2 : |r'-0'| <K,\ v -t'- 0'| > ^\ V-iAl) 



Taking ( = r\ — r and changing the order of integration in £ and r\ by Lemma 16. 1[ we have 

HQ/ 1 (9)ii^. a <cii«iii a / i?(oi 2 / , , ^|m 2/3 - 4 

Jk? 7 {l)eA ( C ):|C|>M} Kl 

\q{r' + <j)' - r))\ 2 da(T')da(cl)')da(r))d( 

Mv) 

<c\\ q \\iz I \m\ 2 m)d$<c\\ q \\i4qf _, a , 

where the last inequality follows from part (i) of Lemma 16.31 with k = 1 and -Ffe(C) is defined in 
(|6.ip . Let us remark that A{rj) C A\(rj) according to the notation in (|3.30|) . 

□ 

Proof of claim (I3.23b[) . We consider the partition Ii{rj) = Ii, a (v) U Ii,b(v) > where 
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h, b (v) :={(e,r,^) €/!(»?) : |„-0|> Jl} ■ 

On the region Ii ta (rf), |£| > ^jj holds. We may follow the same lines as in the proof for the case 
IIk(v) with k = 1, interchanging the roles of the factors q(rj — — £) and g(r + — 77) . As we 
mentioned before, A{rj) C A\(rj) , according to the notation in l|3.30|) and we may apply part (i) of 
Lemma 16.31 for k = 1 . We have 

\\Q Tl Jq)\\w^ <C||g||| 2 ||g||^_i, 2 . 

On the region Ii^(rj), \rj — t\ > -Sk holds. Hence, the proof of the estimate for the region 11(77) 



400 

is valid here, deducing 



hijw^ ^ C||g||| 3 ||<z||^_i i2 • 

□ 



Proof of claim (|3T23c|) 



Let 77 e M 3 \ {0} . The occurrences of q in the term Qi 2 (q)(rj) interact with each other, so that we 
only may bound by the maximal operator just once. We need to consider an extra splitting carried 
out by taking the set of points {6j : 1 < j < AT} in the unitary sphere § 2 , where AT is large enough, 
to get a covering of the sphere T(r]) with AT spherical cups Jj{rj) centered at ftj of radius Ci 1 77I , for 
a certain constant C\ > to be chosen later (with Af ~ ) . We define for every j 



Then 



where 



3=1 3=1 



Rj,(q)(rj) I j i \q(m(v-r)q(T + <f>-i,)\ 



I 7 ?! 3 Jjj(v) J X*(T)JYr,(<j>) 

x \q(r,-4>-0\da(i)da{<}>)da{T), 



and 



X v (t) := {cb E T( V ) : \<f> - r\ > \ V \/40Q , \ V - <j> - r| > |r?|/100} , (3.43) 
Y n {<f>) := U € rfa) : |r? - - CI > b 7 |/100> . (3.44) 

We fix j 6 {1, ... , Af}. In this part, we take an orthonormal reference of M 3 {ei, e2, 63} such that 
ei = 9j, according to the notation used in (|3.42[) . On the integral expression i?^ (g) (77) , we apply 
Cauchy-Schwartz inequality as in (|3.26p - (|3.27[) . For each j, 77 fixed, by Fubini's theorem, we have 

/ / \q(0<l(r + ^-v)\ 2 da(Od<J^)d<j(T) 

J J (ri)Jx v (r)JY v (<j,) 

<[ M)\ 2 I I \q(r + <P- V )\ 2 da(^da(T)da(0. 
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Moreover, 

[ [ [ |9(ry-r')g(^0'-ni 2 ^(e')^(0')^(T') 
Jjjin) Jx v (t<) Jy v (<i>>) 

< [ [ [ ^-r')^-^' -e)\ 2 da{Oda^)da{r') 
Jr( v ) Jr(-q) Jy v (4>') 

\q( V -T')\ 2 d<j(T') f f \q( V -<t>'-0\ 2 M?)d<r(<f>') 
r( v ) Jr(n) Jy v (<i>>) 

<C\\q\\h [ ( \q(T]-cl>'-e)\ 2 da(e)da(<P'), 
Jr( v ) Jy v (4>>) 



bounding by the maximal operator using Lemma 16.21 in the last inequality. Changing the order of 
integration in £ and 77 by Lemma 16. 11 we may write 

\\RjA<l)\\ 2 w e,<C\\q\\ 2 L2 [ \v\ 2 ^ 6 I \q(ti)\ 2 [ I \q(T + <f>- V )\ 2 da(<t>)dcT(T) 

\q(v - 4>' - H')\ 2 d<T(Od<T(<t>')d<T(0dri 

r(rj) JY V (4>') 



A(€) J.l 3 (ri)JX n (T) 



V( V ) JYnW) 



We fix £ G K 3 \ {0} and denote 



Gm--= \vr- b / \q{T + <t>-r,)Yd*{4>)da{T) 

'A(S) J J An) Jx v (t) 

I \q(r,-<t>'-e)\ 2 d<T(e)da(ct>')da(r,). 

r( n ) Jy v (4>>) 

We write r\ G A(£) in cylindric coordinates 77 = £ + sz, with s > and z G {C}" 1- , \ z \ = 1- It is true 
that da(rj) = sdsda(z). Let h(s) := \t)\ = (\(\ 2 + s 2 )^ . We have 

Gj(0= [ I™ Hs) 2f> - 5 I I \q{T + <t>-{Z + S z))\ 2 da{<p)da{T) 

JS 1 JO J Jj-(£+sz) JX €+sz (t) 

x [ [ \q(£ + sz-<f>' - £)\ 2 da(£')da(<j)')sdsda(z) 

<C\\q\\h I r h ^) W - A (3-45) 



S 1 Jo 

\q(Z + sz-(f>' -t')\ 2 do{Oda{4>')sdsda(z), (3.46) 
where the last inequality follows from the following 

Claim 3.1. Let 1 < j < TV, £ G K 3 , « G {C} 1 -, \z\ = 1 tmd s > 0. ffence 



/ / 



\q{r + ^ - + sz))\ 2 da(0)da(T) (3.47) 

) 



L- • 



Proof of claim 13. li 
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We write r, <p in spherical coordinates with respect to the reference {ei, 62,63}: 
£ + sz his) 

H — — (sin^cos^ei +sin^isin(5e2 + cos^ea) , (3.48) 



2 2 
£ + sz h(s) 



(sin£cos7ei + sin £ sin + cos£e 3 ), (3.49) 



2 2 

where tpX <E [0, 7r], S, 7 <E (— n, n]. It holds 

da(<p)da(r) = /i(s) 4 sin ^ sin £ t^y d£ d<$ 

Notice that if r € Jj(£ + sz) then r belongs to the "curvilinear square" from the sphere T(ry) which 
contains the spherical cup </?(£ + sz) given by 

(ip,8)e ^~e ,^+e a x[-£ ,£o], 

where Eq — £o(Ci) satisfies sineo = 2C\. For each £ £ [0, 7r] and ^, we define 

X*(C,1>,S) :={ 7 e (-tt.tt] : ^e%-M} 

7 £ (— 7r,7r] : — (1 — — ■ ) < sin ib cos S sin £ cos 7 

5000 

+ sin ip sin <5 sin £ sin 7 + cos ^ cos £ < 1 



80000 

The integral expression (13.47)1 is bounded by 

C / / / / l ( S ) 4 sin^sin£|g(A(j, S ,^,5,£,7))| 2 d7rfC^^, 

Jf-eo J-eoJO Jx*(C,i/>,6) 

where 

A(j, s, ipi 8, C) 7) := ((sin 1/1 cos <5 + sin£ cos 7) ei 

+ (sini/> sin S + sin£ sin 7) e2 + (cos?/' + cos£)e3). 

For technical reasons we divide the domain which corresponds to the angles ip,S, £,7 into two 
pieces A\ , A2 ■ 

A 1 :={(?/>, <5,£, 7) : |cos(^-£)| < 1 - lO" 9 } , 

A 2 := {ftM,C,7) : |cos(^-£)| >1-10- 9 }, 

where (tp, 5, £) G [ f — £0 , f + £0] x [ — £o ; £0] x [0, 7r], and for each tp, 5, £ fixed 7 belongs to the set 
X*((, tp, S). In this way (|3.47[) becomes bounded by 

C {IIIL + III ) sua V sin C l$(AC/, *, C» T))| 2 <*T dC . 

Remark. We choose A/" large enough in order to take the radius of the spherical cup Jj (77) with 

El 
2 • 



Ci < 10 



Estimate for the domain A\ . 

If Ci < then I cos?/>| < 2C X < 10~ 5 . We have 

1 — 10~ 9 > I cos(t/> — £)| > | sin^ sin£ | — | cos ^ cos£ | 



> \J\ - 10- 10 sin£ - 10" 5 |cos£| 



2(1 
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V 7 ! - 10- 10 y/1 -cos 2 C - 10" 5 | cosC | , 



and 



cos 2 C + 2 • 1(T 5 (1 - 1(T 9 )| cosCI + (1 - 1(T 9 ) 2 - (1 - 10" 10 ) > , 
where the discriminant of the corresponding quadratic polynomial in | cos £ | is 

A = 4 • 1CT 9 (1 - 10" 10 )(2 - 1CT 9 ) > , 

in such a way that £ satisfies 



-5 



| cos C | > -(l - fir 9 )icr 5 + v/10- 9 (1- 10- 10 )(2- 10- 9 ) > 3 • fO 

Keeping in mind that | cosip\ < 10~ 5 , | cos£ | > 3 • 10 -5 , it is true that 
| sin C cos ip cos(S — 7) — sin ip cos £| 

> I sin?/> cos C| — I sin£ cosip cos(<5 — 7)| > | sin-0 cosCI — | sin£ cos^l 



> 3 • 10" Vl - 10^ 10 - 10"Vl-9-10- 10 - 2 • 10~ 5 . 



(3.50) 
(3.51) 
(3.52) 



For each j, s, S fixed we take the change of variables (£, 7, ip) — > A = (Ai, A2, A3) , given by 



A = ((sini/) cos S + sin£ cos7)ei 

+ (sini/; sin 5 + sin£ sin7)e2 + (cos-0 + cos^ea) 



(3.53) 
(3.54) 



We have 



<9(Ai , A 2 , A 3 ) 



h(s) 3 sin C I sin £ cos ip cos(S — 7) — sin ip cos £| 



By Fubini's and Toneli's theorems and estimates (|3.50[) . (|3.51l) . (l3.52p . we may write 

/i(s) 4 sin ip sin ( \q{A(j, s,ip,SX,l))\ 2 d'y d( dS dip 



A, 



< C 



h(s)\q(X)\ 2 d\d6 



Ch{s) / |g(A)| 2 dA. 



Estimate for the domain A2 ■ 



Now we apply the change of variables (ip , S , 7) — > A = (Ai, A 2 , A3) , given by (13. 53|) - (13. 54p for 
each j, s, £ fixed. It holds 



d(Ai , A 2 , A 3 ) 



d(tp, 5, 7) 



/i(s) sin sin £ |sin(<5 — 7) | 



On one hand, since \ip — ?| < e^, we have that | cos V>| < 2Ci. Hence the sinus of ip is lower bounded 
by a strictly positive constant. 

On the other hand, since 7 € X*(£, ip, S) we know that the expression 

sin ip cos 6 sin £ cos 7 + sin ip sin 5 sin £ sin 7 + cos ip cos £ 
= cos(0 — £) cos(5 — 7) + cos ip cos£ (1 — cos(<5 — 7)) 
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takes values between -(1 - 1/5000) and 1 - 1/80000. We have 

1 ~ > |cosO - C) cos{5 - 7) + cos ip cosC(l - cos(5-7))| 

> I cos(t/> — C) cos(<5 — 7)1 — I cosip cos C| (1 — cos(<5 — 7)) 

> (1 - 10 9 ) I cos{6 -7)| -2C X I cos C| (1 - cos((5 - 7)). 
Provided that cos (S — 7) > it holds 

1_ 80500 " ( 1 - 10_9 + 2C il cos CI)cos( ( 5- 7 )-2C 1 |cosC|, 



and hence, 



, 1 - 1/80000 + 2Ci I cos C| 1 - 1/80000 + 2C*i 
cos(5 - 7) < — q , ^ r— Ti ^ 



l-10- 9 + 2Ci|cosC| " l-lO" 9 
Nevertheless, if cos(5 — 7) < we may write 

1 " 80u00 " ~ (1 ~ 10 ~ 9 " 2Cl1 C0SCI) C ° S{6 ~ 7) " 2Cl1 C0SCI 
> -(1- 10" 9 -2Ci)cos(<5-7)-2Ci|cosC|. 

The choice C\ < allows us to write 

1 - 1/80000 +2Ci I cos C| 1 - 1/80000 +2Ci 
C0S{S ~ 7) ^ 1 - IO- 9 - 2d " " l-10-»-2C 1 ' 

and 

vi 1 1/80000 + 2Ci 

cosld — 7) < - — — . 

1 V nl ~ l-10- 9 -2Ci 

Hence | sin((5 — 7)! is bounded below by a strictly positive constant which only depends on C\. 
By Fubini's and Toneli's theorems, we have 

/i(s) 4 sin ^ sin C \q(A(j, s, tp, 6, C, 7))| 2 G?7 d( dS dip 

<C I I h(s)\q(\)\ 2 d\d( 
Jo Jr 3 

= Ch(s) [ \q(X)\ 2 d\ 7 
and claim [3~T1 follows. 

□ 

We return to the expression (|3.45p - (l3.46p . In (|3.46p we write the variables £',</>' in spherical 
coordinates as we did in (|3.48[) . (|3.49[) : 

C + sz h(s) 

C' = - — — 1 — (sin 9 cos 9 a + sin sin 9 e-i + cos 9 63) , 

£ + sz h(s) , . ^ , . j., . 1 j . 

<p = — 1 — (sin ( COS7 ei+smC sm 7 e2 + cosC 63), 

where 9,C' G [0,tt], 0,7' G (-tt.tt]. It holds 

da(C)da(ct)') = h(s) 4 sin 9 sinC' dd dO drf d£. 

For each 9 G [0,tt] and C', 7', we define 

Y*(Q, C, i) := {9 G (-7T, tt] : C' G Y i+sz (4>')} 
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{# G (— 7r,7r] : sin 9 cos 9 sin g cos 7' 
+ sin 6 sin 9 sin f sin 7' + cos 8 cos g > - (l - 1 /5000) } . 



(3.55) 
(3.56) 



(j3~45>(j3~46l) is bounded by 



C\\q\ 



L- 



S 1 Jo 



h(s) 



2/i 



7T />7T />7r 



sin 9 sin g 



it Jo Jy*(e,c,7') 



x |g(B(j, s, 6, 0, C', 7')) \ 2 d9 dO drf dg s ds da(z) 



C\\q\ 



L- 



h(s) 



7T {*7T p7T 



sin 9 sin g 



/o J-ttJq Jy*(e,c,7') 
x \q(B(j, s, 9, 9, g, y))\ 2 d9 dO d-y 1 dg s ds , 



where 



B(j, s ,e,9,c,i') ■= 



-h(s) 



((sin 9 cos 9 + sing cos7')ei 



+ (sin9 sin# + sing sin7')e2 + (cos9 + cos£') e 3)- 
Next for each j, g , 7' fixed we change variables (s, 9,0) — > fi = (yUi, /ii2, ^3) given by 

H = B(j,s,Q,6,g,l'). 
The Jacobean of this transformation is given by 



<9(Ml,M2,M3) 



8(8,9,0) 



s h(s) sin 9 



(1 + sin9 sin£' cos 9 cos 7' 



sin 9 sin£' sin 9 sin 7' + cos 9 cos£') . 



Notice that the mentioned change involves an expression for s in terms of [i which depends on the 
parameters j, g, 7'. Hence, the function h has the same parametric dependence 

h(s) = h(n,j,g,Y). 

Nevertheless, for 9 e Y*(9,g,~/') it holds 

h(s)~\n\ = \B(j,s,9,0,g, 7 ')\ 
thanks to the condition stated in (|3.55[) - (|3.56[) , where 

IM*)I 



(3.57) 



\B(j,s,9,0,g, 7 ')\ 



(1 + sin 9 sin g cos 9 cos 7' 
sin 9 sing sin 9 sin 7' + cos 9 cosg)' 



Indeed, property (|3.57[) is a reminiscence of the condition j^L < | 7 y_^'_^'| < |^| held for g S ((/)'), 
see (pT44|) . 

We conclude that 



Gj(0<C\\q\\ 2 L i 



tf^qi^dndi dg 



10 J-TT. 

= c\\ q \\h I l/f^ 1 !^)! 2 ^, 



and then, we have proved the estimate 



\\RjM)\\w^^<ch\\^- h Aq\\U\q\\^ 
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and since J\f is an universal constant we also have 



\\QiM\\w^ < c k\\ w -hM\l4q\\^-i, 

This ends the proof of claim (|3.23cP and estimate (I3.23[) . 



□ 



Proof of estimate (|3.25|) . This case is inspired on the method used to control the piece Q'jj(q) 
of the cubic term from the Neumann-Born series in the three-dimensional case in |28] , 

Let us start by decomposing the set IV (rj) as follows: IV (rj) C IV<{rj) U IV>(rj) , where 

iv<(v) ~U^,4>)€iv(v) ■■ \tl\r\M< (^ + 4=1 , 



50 V2 

IV>{ri) :=i(t,T,cl>)eIV(r)) : \ V - f|, \q - r|, \r, - <f>\ < ( 1- + ± ) \ v \ 



50 V2 



In fact, if |£| < \r]\ hence 



and if |f | > 4- \r)\ then |r? - £| < A- \r}\ , and it holds 



\V ~ 4>\ < \V - £l + If - <t>\ < (4= + M , 1*7 - r| < |r? - 0| + |«£ - r| < (J^ + -J= 

Taking the changes of variables £ = r\ — t = 7] — t', = r; — 0' in the integral 

|9(0«(»7 - " 4>')q{<t>' ~ e)\da{Od<T(T')da(<p') , 



we notice that Qiv > (q)(v) = Qiv < (q) (t]) , and then Qiv(q)(v) < 2 Qiv < (q)(v)- We take another 

+oo 

decomposition: JV<(?7) C [J (lVk(r)) U IVk(vfj 7 being 



fe=i 

J%fa) == {(e,r,0) G JV<fa) : |0| < |f| ~ 2" fc |r,|} , 
JVfcfo) :={(e,r,0)G/V<(»;) : |f| < \<f>\ ~ 2- fc M} . 



It holds 



(?)IU..) > 
fe=i 

hence estimate (|3.25p follows from the following claims: 

WQivMWw^ < c ^ £k \\q\\ w -\Aq\\w-^Aq\\LA\q\\ w ^i +e ^ , (3.25a) 
\\Qiy k (q)\\w^ < C2- £k \\q\\^_ i: M\w^M\q\\L4q\\ W8 -i + s, 2 (3.25b) 

+ C2^ k ||?||^_ 4 _., 3 ||«||i a |M|^_ i+ .. a , (3.25c) 

provided that e > . 



24 
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Proof of claim (|3.25a)l . By Cauchy- Schwartz inequality, 

1 

w 



QivM{ri)<^ [III \qm(v~T)\ 2 da(Oda{T)da{4>) 

iv k (n) 



\q{T> - c/>')q(<j>' - e)\ 2 d*(e)da(T')M<l>') , 
IVuiv) ) 

where 

IVk(v) = {(£, r, 0) e T( v ) 3 : \cf> - f | < |r/|/100, \4>-t\< It/1/100, 
laM.M < (1/V2 + 1/50)17,1, H < |C| ~ 2- fc H}. 
For each r\ € R 3 , € r(??) fixed, using as above the maximal operator we have 

\q(r'~cb')\ 2 da(T')<C\\ q \\l 2 . 

r(»f) 

Since a(T(r))) = n \rj\ 2 we also get 

\\QivM\\ 2 w^<C\\q\\ 2 L . I H 2/? - 4 / l?(0| 2 / |g(r/-r)| 2 da(r)da(e) (3.58) 

l<^'-o 2 M£')<M0')*7 ( 3 - 59 ) 

c k (n) 

<C\\q\\ 2 ^ \\q\\\ 2 j \ V \ 2 ^ 4 [ \q(0\ 2 (3.60) 

x2- k \ V \ 2 f \q( V -T)\ 2 da(T)da(0dv (3.61) 



2/3-2 



<C2- 2fee || g ||^ , 2 || g ||i 2 / |g(0| 2 / H 



1^1 



2? 



Kl 



2£ 



where 



(3.62) 



/ |gfa-r)| a dtr(7-)dcrfaK (3.63) 



<C2- 2te ||q||^ i , 2 ||q|| 2 / _ e , 2 ||g||| 2 || 9 ||^_, +e , 2 (3.64) 

:= € r(r/) : |£| ~ 2- fc | ?7 | , |f| <(L + ±^\ |„| 
B«fa) == e r(r/) : |£ - r| < M , | T | < (1 + _L^) |,|| , ( 3 .65) 
Ck(v) ■= Ue, 0') € T(, y ) 2 : |0' - e'| < J^L , |^| < |£'| „ 2 - fc h| 

A*(0 := [v e A(0 : |£| < f 1 + |r?|| . (3.66) 

The estimate (EH) allows us to estimate (13381) - (13351) by (|3~BD|) - (i3~M|) . The step from (|3~BuT) - (|3~HT|) 
to (|3.62p ~ i|3.63|) follows from the property |£| ~ 2~ fc |?7| and the change of the order of integration in 
£ and r\ through Lemma I6TT1 Finally, we bound (|3.62[) - (|3.63p by (|3.64[) applying part (it) of Lemma 
IO 

□ 

Proof of claim (|3.25b[) - ([3T25c|) . 
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We split IV k (r]) = IV Ka (rj) U IV k , b {rj) , where 

IVkAv) ■= {(Z,T,cf>) S fV k {ri) : |£| > 2 fc 2 1 ^ | } 

= {((, r, 0) S T(ry) 3 : |0 - £| < h|/100, |0 - r| < |tj|/100, 

|r|, |0| < (1/V2 + 1/50)M, 2- fc - 2 M < |e| < H ~ 2" fc M} , 
^uM := {(£,r,0) G ZV fc fa) : |£| < 2- fe - 2 |r;|} 

= {(^T,0)er(77) 3 : |0-£| < h|/100, |0-r| < |t?|/100, 

|£|, |r|, |0| < (1/V2 + 1/50)M, |C| < |0| ~ 2- k \ V \, |£| < 2- fc - 2 | ?? |}. 

On the domain IVk, a (v)i 2~ h ~ 2 \ r l\ < l£l < 101 < 2~ fe+1 |??| holds; hence following the steps of the 
proof for the domain IVk(r)) we arrive at 

WQfv h ,M\w^ < C2-^ \\q\\ w ^, 2 \\q\\ W -^ 2 h\\L4q\\ w ,-i + s,2 ■ (3.67) 
On the domain IVk,b(v) we have |£ — 0| > 2 fc 2 1 ^7 1 , in fact 

le - 01 > 101 - lei > 2- fc - 1 hi - 2- fe - 2 M = 2- fe - 2 hi . 

In this case, we can bound ||Q/y fc (9)11^^ 2 by a similar expression to (|3.58[) - ()3.59|) replacing 
B k (r)) by the set {£ € r(7?) : |£| < 2- fc - 2 M, |f | < (i + ^=)|r?|} and the domain C fe (rj) by the set 

(£', 0') € r(r,) 2 : 2- fc - 2 h| < |0' - e\ < |^ , |£'| < |0'| ~ 2- fc h| 

By the estimate (|3.35[) , changing the order of integration in £ and r\ by Lemma 16. 1[ multiplying 
and dividing by j^p^ and applying that < 2 2e (~ fc_2 ) , and finally by part (ii) of Lemma 15751 we 
obtain 

nr '<iv k jb wnw.2 iiyii v ^-i- e ,2iiyiiL2||y|i vt/3 _ : L 
The expressions (|3~67f and (f3768]l lead up to claim (|3.25bj) - (j3.25c|) 
This ends the proof of estimate (|3 . 25[) and Proposition [3j 



Qfv. JMwe* <C^ k \\q\\ 1 e , 2 ||g||| 2 ||g|| .„ i +e>2 . (3.68) 



□ 



4. Proof of Theorem \S\ (remainder term a < n/2). 

The control of the remainder term R; , where I is as in the statement of Theorem [SJ follows from 
the next proposition by choosing Cq large enough in (jTTTJ) . For Co > (2||<7|| ^0,2)"* we obtain the 
convergence of the series 



+00 

E 

3=1 



, 1 

2C 4 Nlw° 



Proposition 4. Let 71 € {2,3}, q £ W a ' 2 (W n ) compactly supported and < a < n/2. Assume that 
Co > 1 7 j > 4 if n — 2 and j > 5 if n = 3. Then, for any (3 € K suc/i i/iat /3 < a + 1, holds: 

WQMWw*.* < c(a,p)c 1 7/i (2C- i \\ q \\ w ^y \\q\\ L 4q\\w^- t 4 - 1 ) 



:2l i 
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Proof of Proposition [H We write Rg, k f(x) = e^ ke - x R k (e ike <^ f(-))(x). It holds 

QM)®= f e* ke -y(qR k y- 1 (q(-y ke ^)(y)dy= f e 2 ^ '{qRe^ 1 (q)(y) dy , 
where k — |£|/2 and 8 = — £/|£|. Let 7 G R be such that 7 < f3j, being 



Pi ■■= { 



!(i-2) + f(i-l), if a< i and n = 2, 

(i-3)(| + f) + l, if \<a<l and n = 2, 

+ (j - l)f - I , if < a < f and n = 3 , 

'1 1 a\ , 1 ;t 3 x „, / 3 



(j - 3) (| + f ) + I , if I < a < § and n = 3 . 
Taking the change of variables £ = — 2k8 with > , 6 G S*™ -1 , we have 

\\qm\\%^ = / i^i 27 i5^)(oi 2 rfe= / \e 7 x*moM(o\ 2 dc 



r + 00 r- 

<C„2 2 W k 2 ^" 1 / IKgiJe,^'- 1 ^)!!!,^^)*. (4.2) 

By Cauchy-Schwartz inequality, || (qRg^)^ 1 (q)\\ L i < C\\q\\ L 2\\R e ^{qRg tk ) J ~ 2 {q)\\ L 2. Using the es- 
timate given by Lemma 3.4 in j26j for the operator Rg -k and the next inequality for products of 
Sobolev spaces due to Zolesio 

||/g||iv=3-p < C \\f\\w a i-pi \\gWw2.p2 , 

where ai, 03 > 0, < otj, pj > p, j — 1, 2, 011 + a-2 — 0^3 > n(l/pi + I/P2 — 1/p) > 0, we arrive 
at 

\\ReAlRe, k y- 2 {q)\\L* < Ck^\\q\\^l 2 , 
where 7,- := -(j - 1) + ^=±(j - 3) (1/2 - a/n) + 2=1 max{0, | - ^} . All this leads us up to 

where the integral converges if 27 + 27^ + n < 0, that is to say, if 7 < j3j. Notice that Pj = — 7- — jj. 
In this way, we have proved that, for 7 < fy, it holds 

WQMWw^ < (t) 7 ^ ^ L2 ^w^ ■ (4-3) 

Let e = s(a, @) := (a+ 1) - /3 > . Keeping in mind 2/3 = 2(/3,- - e) + 2 (a + 1 - /%) and a + 1 < 
for our j, we write 

IIQ;(g)IU.= <c Q+1 - ft WQMWw*-.' 



\\Q j (q)\\%, t2 <C2 2 ~< fc^-^ifcllgll^lkll^ 



= C(a^)2*C ^||g|U 2 |klfc' 

where the last inequality follows from formula (|4. 3[) in the case 7 = /3j — e . In our setting, 2 l3j < 2- J . 
Moreover, j3 — < a + 1 — y8j < —\j + lp , for our j, and C > 1. We have proved (|4.1[) . 

□ 
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5. The case a > n/2. 

In this section we are going to extend Theorem |H Theorem [5] and estimates (|2.5I) . (|2.6[) . (|2.7|) and 
(|3.ip - (|3.2p for any a > 0. Then Theorem Q] will follow from these estimates for any a > 0. We start 
with a Leibniz' type formula for derivatives of Qj(q) which we state as follows 

Theorem 6. Assume that a g N" , j S Z , j > 2 and Zet g S W^' Q '' 2 (R n ) 6e a compactly supported 
function. Then 

Pi,...,0j>Q 

Remark. From the proof of Theorem [6J one also deduces the formula 

D a QM= E ~al — ~~ — ~r~\ QA D01 1 >■■■■> D^q) , (5.1) 

/3 1+ ...+ft=a^-' 'Pi' 
0i,...0j>O 

for the same hypotheses on q. 



Proof of Theorem [6J Writing the resolvent Rk as the convolution operator with the outgoing 
fundamental solution to the Helmholtz equation (see [BJ, [2T] ) 

^ M -q.^/» £|gggW) , (5.2) 

where C», = ,L ,w, and if, 1 ' „% /n denotes the Hankel function of the first kind and order 
(n — 2)/2 (see [35]), we have 

Q~W)(-2k6) = J e^yqmRmy-^e^m dy 



/ e lfc9 ' '^(an) TT (0fc(«i " xi +1 )q{x l+ i) ) e M ' x * dx . 



1=1 
and 

F(Q j (f 1 ,...,f j ))(-2k6) (5.3) 

, i-i 
= / e ik0 -^h{ Xl ) TT (&(x, -x I+ i)/, +1 (a!, +1 )) e ike - x * dx , (5.4) 

where dx = dx\ ■ . . . ■ dxj and € M n , for any I = 1, j . We know that J- (D a Qj(q)) (—2k6) = 
{-i2k6) a cf~(g){-2k6) . Taking the change xi=x x +y u 2<l< j, it holds 

QM(-2k9) 

3-1 



f e ake - *iq{ Xl ) TT Ufa ~ x l+1 )q{x l+1 )e- lke - ta-'+O ) 

./(K™) J V ' 

ft i~ x 

\ \ e i2k6 -^q( Xl ) TT (Myi-m+iMxi + yi+i)e~ zke ' iv '~ Vl+l) ) dyd Xl 



1=1 

where y\ — and dy = dy^ • . . . • dy?.. Integrating by parts, we have 
(-i2kd) a Q~(q~){-2k8) 



2N 
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(-1) H / / ._ {D^e aM ^)q(xi) J ]\Uk(yi-yi + i)q(xi + yi + i)e- ike <v>-^A dyd Xl 

dydxi 



In 



E 



o! 



/3i+...+/3j=a 
X 

1=1 



0i 



! . 



•ft! 



9(*i) J] (^ fc ( W - 2/i+iMzi + y/+i)e-' ifc0 - to - y,+l) ) 
z=i 



JJ (&( w - yi+x)e- ike - to-vi + i) D Pi + x q ( Xl + w+1 )j dl/daJl , 
where we have applied Leibniz' formula: 



D' 



u •...•/*)= E 



Px,-,Pk>Q 



'V- ... - A! 



/>' ./! • ... • 



Finally, 



Pi+...+P j= a 
Pi,...,P s >0 



e ik6 - Xl D^q(xi) 



II - ^+i)£> ft+1 9(^+i) ) e lfee ^ <fo, 



and remembering the expression (|5.3p - (|5.4|) , we end the proof of Theorem [6] 



□ 



To prove Theorem [T] in case a > n/2 we use induction on [a]. We need to use the boundedness of 
the j-multiple scattering operators, see the notation, when acting on (gi, ...qj) where j — 1 of them 
are equal to q. Namely 

Proposition5. Let n £ {2, 3}, a £ R with < a < n/2 and let us suppose that qi,q2 £ W a ' 2 (R") 
are compactly supported functions. Then Q2(<7i,Q2) G W@ ' 2 (M' 1 ) + C°°(R n ) , for any £ R smc/i 
£/iat < /3 < a + | . Moreover, there exists a constant C(a, 0, qi, q%) > which just depends of a, 
and the supports of qi , q 2 such that 

II Q2 (gi,92)||^,2 < C{a,0,q 1 ,q 2 ) max{||gi||^ , 2 , Ikall^.a} ■ 
Proposition [5] follows by polarization of estimate (|2.5p . 

Next propositions [6j [7] and [8] follow from the proofs of the analogous estimates (|2.6p . (|2.7p . (I3.ip ~ 
(|3~2")1 and Proposition H 

Proposition 6. Let 71 € {2,3}, a £ R wi£/i < a < n/2 and 51,(72,93 as gi /rom Proposition^ 
Then Q 3 (gi,«2,?3) € H" 3 ' 2 ^") +C°°(R™) 7 /or any /3 € R ftoZdm^ < /3 < a + 1 i/ n = 2 and 
< < a + 1/2 if n = 3. Moreover, there exists a constant C{a, (3, gi, 02, 53) that just depends of 
a,0 and suppgi, suppg2, suppg 3 sttc/i i/ia£ 



IIQ3(gi,92,93)||^,: 



(5.5) 
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<C(a,{3,q 1 ,q 2 ,q 3 )(^2 ll? CT (i)IMka(2)|Ml ff a(3)||^c«,2 (5.6) 

+ 11^(1)11^-1,211^(2)11^-^211^(3)11^,2 (5-7) 

reS 3 

+ lkc;(l)ll v ^-^- e ,2||Q'a;(2)||z, 2 lka;(3)ll^.0i ( 5 - 8 ) 

where e — a + 1 — (3 > if n — 2 and e = a + i — (3 > if n = 3 . 

Proposition 7. Let a £ R with < a < 3/2 and qi, 02, 03, 04 as gi /rom Proposition^ for n — 3. 
ThenQ A (q 1 ,q 2 ,q 3 ,q4) G W^' 2 (R 3 ) + C°°(R 3 ) , /or any /? G R withO< [3 < a + 1/2 . Moreover, i/iere 
exists a constant C(a, /?, gi, 52, g3, 04) > jusi depending of a, [3 and the supports of q±, q 2 , 03, 04 suc/i 
£/ia£ 

||Q4(gi) 92, 93, 04)Hvi/,3,2 

< C(a,/3, gi,g2,g3,04) ( ||9,T(l)||L2||g (T( 2)||L2||g (T( 3 ) || i 2||g (T( 4 ) || l y Q ,2 (5.9) 

+ lkr(X)ll^-J,2||«r(2)|U»lkT(3)IU a lkr(4)|lvfcro,,a (5.10) 

t<ES 4 

+ X! 11^-^,2 lkcc,(2)ll^- e ,2|ka;(3)IU 2 lka,(4)ll^ a ,2 (5.11) 

+ lkp(i)ll^-i- s ,2||flp(2)IUHkp(3)IU 2 ll o p(4)llw«.0) ( 5 - 12 ) 

where e = a + j — (3 > 0. 

Proposition 8. Let us assume that n £ {2, 3} , a <E R, < a < n/2 , q\,...,Qj £ W a ' 2 (M. n ) are 
compactly supported functions and Co > 1 . Hence, for any (3 G R suc/i that (3 < a + 1 : 

||Qi(«l,...,«,0IUw < C(a,/3)C 17/4 (2Co* max IMI^) , (5.13) 
where j > 4 if n = 2 and j > 5 if n = 3. 

6. Appendix. 

In this section we state two results, Lemma 16.11 and 16.21 which are often used in this work and 
state and prove an important result, Lemma [6731 in order to demonstrate Proposition [3l 

Let V be the submanifold of R 2n V := {{n, f) G R™ x R™ : £ • (£ - rj) = 0}. Then V can be viewed 
as a bundle of spherical sections V = {(n, £) G R" x R™ : 77 G R™, £ S r(n)}, or as a bundle of plane 
sections: V = {(r),£) £ R™ x R" : £ G R" , 77 G A(£)} , where L(n) and A(£) are dehned in Q3J) and 
(jl.5p . In this context, the following lemma from allows us to change the order of integration in 
£ and r/ . 

Lemma 6.1. Let V = {(?/,£) G R" x R™ :£•(£- 77) = 0}. Let dcr r) (£) oe tfie measure on L(n) 
induced by the n-dimensional Lebesgue measure d£ and Zei da^(rj) be the measure on A(£) induced 
by the n-dimensional Lebesgue measure drj. Then 

In] 

da v {Qdr) = j-r der c (n)d£. 
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The following lemma in [28] is used several times in this work. 
Lemma 6.2. Assume that the support of q is contained in the unit ball. Then we have: 

(1) // e, £' £ R» satisfy \£ - ?\ < 3, then \q(£)\ < CMq{C) . 

(2) \\q\\ L <» <C\\q\\ L 2. 

(3) For < 7 < \\q\\-fr-i,2 < C||g||i2 , where C depends on the size of the support of q. 

The following lemma is fundamental to control the spherical term Q(q) of the quartic term. 
Lemma 6.3. Let (el 3 \ {0}, f3 G K , e > 0, k G N. We denote 

F k (0:= f M 20 -' I [ \q(v-^-a 2 Me)da^')Mv), (6.1) 

G(0:= [ |^- 2 + 2e / -T)| 2 da(r)rfaC7), (6-2) 

Ja*(£) Je^) 

where Ak(rj) , B^(rj), A*(£) are defined in (|3.30p . ()3.65|) . (|3.66p respectively. Then 

(i) Fk(£) < C 2~ 2k L 3 |A| 2 ^ _1 |g(A)| 2 dA , /or some constant C independent of £ , fc and g. 
(u) G(£) < C J R3 |A| 2 ^~ 1+2£ |g(A)| 2 ciA , /or some constant C independent of £ cmd g. 

Proof of Lemma 

• Proof of (i). For rj £ A(£) we write 77 = £ + sz and /i(s) := |?7| = (|£| 2 + s 2 ) 2 , where s > and 
z G {C}^ ? M — 1 ■ F° r simplicity, we don't specify the dependence of variables with respect to £ 
since it is fixed along the proof . It holds da(rj) = sdsda(z) , where der(z) denotes the measure on 
the unitary circumference S 1 in the plane ■ We have 

Fk(0= [°° f h(s? fi - A I I m + sz-tf -a\ 2 d<y^)da( ( f>')da(z)sds. (6.3) 

JO JS 1 J JA k (s,z) 

Fixing z, s we parametrize (jJ G T(£ + sz) by v, u G S 2 , respectively: 

£ = — 5 h_ 2 _l ' — 2 ~2~ U ' u '' yeS ' 

where dcr{(J) — C h(s) 2 da(v) , da(ip') = C h(s) 2 da(u) . The domain of integration for v, u is given 

by 

{(v, u) G § 2 x § 2 : \u-v\< 2~ k+2 ,l + u-v> 1/5000} , 

since < 2- fe+1 /i(s) implies |u-i>| < 2- fe+2 and |£+sz-0'-£'| > ^ implies > ^_ . 

Let {-Dj,/t : j G {1, ... , iVo2 2fe }} be a finite overlapping cover of the sphere S 2 with overlapping 
constant independent of k such that D is an spherical cup of diameter and Nq an appropriate 
constant. For each j we define 

Dj, k := I u G § 2 : \u - v\ < 2~ k+2 ,l+u-v> — \- , for some v G D jik 



The expression (|6.3|) is bounded by 

N 2 2k „ „ f+oo 



c E / / / /. *(- 



\2/3 



2 



da(u) s ds da(v)da(z) . (6-4) 
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Notice that for every j € {1, ... ,N 2 2k }, u G Dj k and v £ Dj :k , we have |w + v\ > ygg. In fact, 
since u € Dj ik there exists a v' G Dj ik such that 1 +u ■ v' > 1/5000 , and hence \u + v'\ > 1/50. We 
have 

\u + v\ > \u + v'\ — \v — v'\ > diamDj.fc 

50 

1 2- fc 1 1 1 
~50~l^-5l)~100~ 100' 
We take spherical coordinates for u with respect to the canonical reference of M 3 , 

u = (cos 9 sin rp , sin 9 sin rp , cos ip) , (6-5) 

with dcr(u) — sintpdipdO . We bound (|6.4p by 



/l(s 



Is 1 Idj^ Jo J 


' L h[s)2P 







sin i/> d^)d9 s ds da(v)da(z) 



where u(if),9) is given by (|6.5[) and ZA* fc := {(^, 9) e [0, 7r] x [0, 2tt) : u(?p,9) e Dj_ k } . We remark 
that for (ip, 9) G D* fc it holds 

\u^,9)+v\>^. (6.6) 

For z G S 1 , j G {1, iV 2 2fe } and u e Dj, k fixed we consider the change (s, 0, — * (Ai, A2, A3) = A 
given by 

A = £ + S z - 0' - e = (uty, 0) + «) 

= — ( cos sin ip + ui , sin sin V 1 + u 2 , cos ^ + «3 ) , 

where u = (i>i, i>2, ^3). We have dX = sfe ( s ) sin V'|i+"W',fl)-t'l dsdipd9. From condition (|6.6p we deduce 
that |A| ~ /i(s) . We define this family of sets with overlapping constant independent of k contained 
in the interior of convex cones in R 3 : 

H jtk := {r(u + v) : u G D jik , v G * , r < 0} . 

Since a(D 3 , k ) ~ 2~ 2fc , we have 

F k (0<cY, / / lAI 20 - 1 !^)! 2 ^^)^) 

f /A r o2 2fc \ 



< C2- 2fc / |A| 2/,_1 |§(A)| 2 dA. 



□ 



• Proof of (ii). We also express 77 € A*(£) as 77 = £ + sz , with s > 0, z G {£} j M = 1 • 

We use the same notation h(s). Since |£| < -f |^| ; it is true that s > 0.9|£| . It holds 

dcr(rj) — sdsd<r(z) . For z £ S 1 fixed, we take a reference {ei, €2, e^} in R 3 such that z — . Then 

£ = (^1)6)0) and r] = (Ci,£,2, s) . We write for r G + sz), 

£ + sz ft,(s) „ 
t = ^— + « , v G S 2 , 
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with da(r) = C h(s) 2 da(v) . We express v in spherical coordinates with respect of our reference: 
v = ( sin 9 cos 8, sin 9 sin 8, cos 9), < 9 < 7r , —tt<8<tt, 

where da(v) = s'mQdQdS. For s,z fixed, let 

M(s, z) := {(9, 8) € [0, tt] x [-7T, ?r) : r(s, z, 9, (9) e^(^ + sz)} . 

We obtain 



G « )s /r // 

JS 1 J0.9\t\ J JA 



s h(s) 



2/3+2e 



sin 9 



IM(s,z) 



q ( — (£i — /i(s) cos 8 sin 9, £2 — /i(s) sin 8 sin 9, s — /i(s) cos 9) 



dQdOdsda(z) 



For each u, we take the change of variables (s, 9, t 



A= (Ai,A 2 ,A 3 ) given by 



1 



A = 77 — r = - (£1 — /i(s) cos 8 sin 9, £2 — h(s) sin 8 sin 9, s — h(s) cos 9) 

It holds d\ = ^ sin 9 |s - /i(s) cos 9| ds dB dd. 

Since r S S^(£ + sz) we know that |r| < f y= + \r]\ and hence, 

1 



1 



1 



y/2 50 



H < \v-r\< \v\- 



Since (9, 8) S .M(s, z) it holds |A| ~ /i(s) ( \r] - r| ~ |??| ). Moreover, |£ - r| < ^ and the angle 7 
between 77 — £ and r] — t satisfies | C0S7I > C > 0. That is, \(rj — £) • (77 — r)| ~ |?y — £| |r? — r|. This 
says that |s — h(s) cos 9| ~ | A| , since z • (rj — r) = s ~ h ( s ) cose _ On our domain of integration we have 
\s — h(s) cos 9| ~ s: 

\s — h(s) cos 9 1 ~ |A| ~ h(s) ~ s, 
where the condition ft,(s) ~ s follows from s > 0, 9|£|. 



We conclude 



G(0<C[ f \q(X)\ 2 \X\ 2 ^ 1+2E dXda(z). 



□ 
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